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We consider the evolution ol a bulk scalar field in anti-de Sitter (AdS) spacetime linearly coupled 
to a scalar field on a de Sitter boundary brane. We present results of a spectral analysis of the 
system, and find that the model can exhibit both bound and continuum resonant modes. We find 
that zero, one, or two bound states may exist, depending upon the masses of the brane and bulk 
fields relative to the Hubble length and the AdS curvature scale and the coupling strength. In all 
cases, we find a critical coupling above which there exists a tachyonic bound state. We show how 
the 5-dimensional spectral results can be interpreted in terms of a 4-dimensional effective theory 
for the bound states. We find excellent agreement between our analytic results and the results of a 
new numerical code developed to model the evolution of bulk fields coupled to degrees of freedom 
on a moving brane. This code can be used to model the behaviour of braneworld cosmological 
perturbations in scenarios for which no analytic results are known. 



I. INTRODUCTION 

Branes offer an alternative to the traditional Kaluza- 
Klein approach to dimensional reduction. Some degrees 
of freedom are restricted to branes, while other fields, 
including gravity, propagate in the higher-dimensional 
bulk spacetime. In this case the extra dimensions may 
be non-compact Brane and bulk fields are usually 
treated as independent degrees of freedom, but in general 
one should include interactions between brane and bulk 
fields. For instance, in order to study inhomogeneous 
perturbations about a brane-world cosmology we need 
to consistently include bulk metric perturbations coupled 
to matter perturbations on the brane @, IS EL IS IS 0> @l • 
Once we include a finite coupling between brane and bulk 
fields it is no longer possible to think in terms of a brane 
field without consistently including the associated bulk 
field configuration, and vice versa. 

In this paper we will consider the case of linear cou- 
pling between a four-dimensional scalar field on a de Sit- 
ter brane and a five-dimensional scalar field in an anti-de 
Sitter (AdS) bulk spacetime. We will develop an ana- 
lytic approach based on the resonant scattering modes 
of the bulk field as has recently been applied to the 
study of quasinormal modes of the gravitational field 
about Randall-Sundrum @ and Einstein-static [IS E| 
branes. We identify both evanescent scattering modes 
(quasinormal modes) of the bulk field, which decay at 
future null infinity, and normalisable bound states. Our 
linearly coupled brane and bulk perturbations about a 
de Sitter brane may be seen as a simplified model of per- 
turbations in the inflaton field driving slow-roll inflation 
on a brane which are linearly coupled to scalar metric 
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perturbations in the bulk [HHSHi- 

Linearly coupled brane and bulk fields on a Minkowski 
brane have been previously studied in a Minkowski bulk 
[IS Ell and in an anti-de Sitter (AdS) bulk 16]. In a 
Minkowksi bulk there is always one and only one bound 
state, and above a critical cou plin g the bound state de- 
scribes a tachyonic instability (la ]. In an AdS bulk one 
finds again a tachyonic bound state above a critical cou- 
pling, but below the critical coupling there is no bound 
state (l6| . On the other hand we know there is a bound 
state for a light bulk field in the case of zero coupling in 
the presence of a de Sitter brane. The spectrum of free 
bulk modes is qualitatively different in the presence of 
a de Sitter brane as only modes above a finite mass gap 
(proportional to the Hubble rate on the brane) can escape 
to null infinity [17|. We will investigate the existence of 
bound states with finite coupling between bulk and brane 
fields in the presence of a de Sitter brane, finding zero, 
one or two bound states depending upon the parameter 
values. 

We also present a numerical code to solve for the evolu- 
tion of the bulk field coupled to the scalar field on a mov- 
ing brane in the bulk spacetime. This is a dev elop ment of 
Seahra's numerical code pj| (see also 0, [H, US HH, H2 ) 
for bulk perturbations with an arbitrary brane trajectory 
but no brane coupling, to include coupled degrees of free- 
dom on the brane. We use this to test our analytic re- 
sults, but we also expect this numerical code to be useful 
to investigate a wider range of models such as metric per- 
turbations in a bulk spacetime which have a non-trivial 
coupling to the brane fields and are thus less amenable 
to analytic study [HESE3. 

The layout of the paper is as follows: In Sec. HH we in- 
troduce our model and the notation used throughout the 
paper. Sec. lIIII discusses the spectral solution of the equa- 
tions of motion and the system's resonant excitations. An 
effective 4-dimensional description of the model's bound 
state sector is developed and compared to the full the- 
ory in Sec. HVl Direct solution of the wave equations in 
the time domain via numeric computation is the topic of 
Sec. El Finally, Sec. IVII is reserved for our conclusions. 
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II. STATEMENT OF THE PROBLEM 

In this section, we describe the model analyzed in this 
paper. We consider a portion of AdS space M. bounded 
by a brane dM.^. The bulk geometry can be usefully 
described by two different coordinate systems 



£ 2 

dsl = -^(-dt 2 + dx 2 + dz 2 ) 



( HI 



\ sinh H y 



{-dr 2 + e 2HT dx 2 + dy 2 ). (1) 



The hrst set of coordinates ("Poincare") are related to 
the second set ("brane normal") by 



Ht 



-Ht 



cosh H y, Hz 



-Ht 



sinh Hy. (2) 



Here H is the Hubble parameter: a positive constant 
with dimensions (length) _1 . 

The trajectory of the brane in the brane normal coor- 
dinates is simply 

y = j/b = H _1 arcsinhi7^, (3) 

which gives the induced line element 

dsl = - dl ~ 2 + e 2HT dx 2 = (H^i-d-q 2 + dx 2 ), (4) 

where we have defined the conformal time rj by 



Htj 



-Ht 



1] < 0. 



(5) 



From ([4]), we see that t is the proper time on the brane. 
We can also write down the trajectory of the brane in 
Poincare coordinates: 



t = t b = 7r/, z = z h 



-rjHi, 



where 



7 



= v 7 ! + H 2 ( 2 = cosh Hy h . 



(6) 



(J) 



We write the normal vector to the brane and the brane 
induced metric as 

n a d a =d y = {-H V )(-H£d t +~fd z ), (8a) 
hab = 9ab - n a n b . (8b) 

Finally, note that our bulk manifold M. corresponds to 
y > y b , which implies that t < and z > by the coor- 
dinate transformation ^ . We assume the Z2 symmetry, 
which implies that the "other" half of the bulk (y ^ y b ) 
is the mirror image of M. . 

We consider two massive test scalar hclds propagating 
on this background geometry. One of the fields (j> exists 
throughout the bulk, while the other q is confined to the 
brane. These are governed by the action 



S=~ I {g^dacfrd^ + m 2 ^ 

IM 

' ' (h ab d a qd b q + [l 2 q 2 + 2[3cj ) q). (9) 
dM h 



The bulk part of the action is twice the usual expression 
due to the Z2 symmetry. The coupling constant f3 has 
dimensions of (length) -3 / 2 and can be taken to be pos- 
itive without loss of generality. The masses of 4> and q 
are m and [i, respectively. Variation with respect to each 
field yields the equations of motion 



<TV a V fc = m 2 
h a ^V a V p q = 



li 2 q 



(10a) 
(10b) 



where h a /s and V Q are the intrinsic metric and covari- 
ant derivative on dM^, respectively, and tfib is the bulk 
field evaluated on the brane. The bulk field satisfies the 
boundary condition 



(n a d a (b) h = Uq. 



(11) 



In what follows, we will almost always work with the 
spatial Fourier transforms of the fields 



3 l, _-ik x 



(r,x,y) = j d'-k 
?(r,x) = J d 3 ke- lk *q k (r). 



c(r,y), (12a) 



(12b) 



The k subscript on (f> and q will be omitted in all the 
formulae below, with the exception of those in Sec. HV| 
where we work with the fields in real space. 



III. SPECTRAL ANALYSIS 

In this section, we seek the solution to the equations of 
motion for our system (|10j) via a spectral decomposition 
of the fields. Our main goal will be to find and classify the 
resonant modes in the problem by looking at the poles of 
the scattering matrix. Throughout this section, we use 
the brane normal (t, y) coordinates. 



A. Canonical wave equations and the scattering 
matrix 

To put the wave equations into canonical form we in- 
troduce the re-scaled fields <p an d q, defined by 

<jy{T,y) = e- 3HT / 2 (sinh Hyf/ 2 4>(T,y), (13a) 
qir) = e- 3HT / 2 q( T ). (13b) 



Then, the wave equations (|10|) reduce to 

+ fc 2 e- 2 ^] 4> = [fig - U(y)} 4>, (14a) 
[d 2 + k 2 e- 2HT ] q = [ftf 2 - /J 2 ] q 0(Hi) 3 /^ h , (14b) 

where we have defined the potential 

iJ 2 (15 + 4m 2 ^ 2 ) 



U(y) 



4(sinh Hy) 2 



(15) 
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The structure of the equations of motion (fl4|) suggest a 
spectral decomposition of the form 



(16a) 
(16b) 



Here, T u (r) is an eigenfunction of the temporal operator 
appearing in (|14[) : 

[d 2 + k 2 e- 2Hr ] T u (r) = -cj 2 H 2 T^t). (17) 

This is satisfied if we take 



J — OG 

q(r) = / du)T u (T)qu. 



T u (r) = 



2H 



Hlo 



2 sinh 5 



1/2 



3-i 



k 

— ( 
H 



-Ht 



(18) 

where J„(x) is the Bessel function of the first kind. We 
note that T u reduces to a single-frequency plane wave in 
the limit of large r, 



!/ 2 e iu)Hr 



2sinh7raj/ r(l — iu>) ' 



[Ht > hxk/H]. 



(19) 

The normalisation coefficient in (TTH)) has been selected 
such that 



/+oo 
drTUr)T:,(T). 
OO 



Putting dT6J) into (JUjl yields 
= [d 2 - U(y) + uj 2 H 2 ] ^ 



(20) 



(21a) 



= [uj 2 H 2 + \H 2 - M 2 ]4, - (3(H£) 3 / 2 4> hw . (21b) 



The general solution of (|21ap is 

My) = a t<t>t(y) + a JC(y): 

^(2/)=AA + (sinhi7y) 1 / 2 ^_ 1/2 (coshi7y), 
C(2/)=^(sinhi/ 2/ ) 1 / 2 g,« ) _ 1/2 (coshi7y), 

where 

Tre*™P?(x) - Q%{x) sin7r(cr + a) secw 



(22a) 
(22b) 
(22c) 



BZ(x) 



M ± = 2e~ i7rQ ( 



r((j + 3/2)r((7 + i/2) 

H smn 7tcj / 



= \/4 + m 2 £ 2 . 



(23) 



Here, P"(x) and <5"(x) are associated Legendre func- 
tions .23J, 8.7]. The coefficients TV* have been selected 
such that when we use the standard asymptotic expan- 
sions [23l, 8.776], we have 

[Ht > In k/H and Hy > 1] . (24) 



Hence, </> reduces to a superposition of incoming and out- 
going plane waves in the "asymptotic region" defined in 
the square brackets; i.e., far into the future and far from 
the brane. 

We have yet to enforce the boundary condition (jlip . 
In terms of the mode decomposition, this reads 



y(Ht) 



-3/2- 



(25) 



When $22]) is substituted into (|21b|) and ([25]), we ob- 
tain a homogeneous system of two equations linear in 
(ajjjQw) which can be solved for the amplitude ratios 
a t> The derivatives of Legendre functions appearing 
in this system can be simplified using recurrence relations 
[H, 8.732]. 

We define the scattering matrix (i.e., reflection coeffi- 
cient) Suj as the ratio of outgoing to ingoing flux in the 
asymptotic region (cf. Eq. I24[) at a particular frequency 
uj. Explicitly, we obtain 



c _ a " 
— ~rp 



7V+ 



wT 2 (iuj)N- 



cos7r(a + iui) 
sinh TTuj 



+ iire 



c+^ +1/2 (7) + c-^_ 1/2 (7) 



c +QL + i/ 2 (7) + C-QL-i/ 2 (7)_ 



where 



c+ = 



9 
4 



EL 

H 2 



a 



iuj 

2 



ip 2 
2H 2 



, 2 y M 

+ «7 H At 

4 if 2 



(26) 

(27a) 
(27b) 



Before moving on, we should mention cases for which 
(|26[) is not a valid expression for the scattering matrix: 



iuj - 1/2 ± a ^ Z, 



Z, 



(28) 



where Z represents the integers. The first of these re- 
strictions is in place to ensure that the mode functions 
(f)^ are linearly independent [23|, 8.707], while the second 
is required for the asymptotic expressions (fM]) to be ap- 
plicable [23|, 8.776]. When to is real, both conditions are 
automatically satisfied. 



B. Resonances 

Intuitively, a resonance of a given system is a mode of 
oscillation where the response ("output") is maximised 
relative to the stimulus ( "input" ) . In scattering theory, 
this implies that the outgoing flux associated with a res- 
onant mode completely dominates the ingoing flux. This 
leads to the formal definition of a resonant excitation: 
A mode whose frequency corresponds to a pole of the 
scattering matrix S u when u> is continued into the com- 
plex plane [111, HH HH, [2(| . The physical importance of 
these resonant modes is that they describe the dominant 
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behaviour of any given system, regardless of the choice 
of initial conditions. In this context, we expect that the 
late time behaviour of the brane and bulk fields to be 
well approximated by resonant solutions irrespective of 
the configuration on some initial data hypersurface. In 
this subsection, we study the poles of in some detail. 



1. General properties of quasinormal modes 

When the scattering matrix (|26p is continued into the 
complex plane, one discovers that it has a number of poles 
on the imaginary ui axis that violate the restrictions (j28|) . 
Since the scattering matrix itself is undefined at these 
locations, these poles are ignorable and we omit them 
from the subsequent discussion. We can show that all 
of the physical poles that remain correspond to roots of 
1Z U , where 



t 1 2 9 



i 

a — iuj 

2 



7 I - - 



liO 

2 1 



1(7) 
2 



2 9\ Iff 
H 2 U 4/ 2H 2 



Q a 1(7). 



(29) 



This condition that TZ U = is the same as demanding 
that a particular mode function satisfies the boundary 
condition at the brane and is a purely outgoing wave in 
the asymptotic future; i.e., a+ = 0. Hence, the resonant 
states we have defined here are akin to the quasinormal 
modes of black hole perturbation theory [2(| . In Fig. [TJ 
we show a plot of \R.~ l \ as a function of complex fre- 
quency for a certain choice of parameters. 

It is useful to note that the full spatial and temporal 
dependance of a resonant solution for the canonical and 
physical bulk fields are 

4%(t,v)=T u (t)$ u {v), (30a) 
05(r, V ) = e- 3HT / 2 ( S mhHy)V 2 T u (T)My), (30b) 

respectively. In both cases, ui is a root of Tt u . The asymp- 
totic behaviour of these resonant solutions is crucial to 
the discussion below. We have for late times 



lim 

Ht — >oo 

lim 

Ht — >oo 



%(T,y)\oce- rH \ (31a) 
^^lloce-^*; (31b) 



while far from the brane, we obtain 

£(r,y)|oce ri ^, 



lim 

Hy — *-co 



lim 

Hy — >oo 



(31c) 
(31d) 



1 3 

£ -10 




FIG. 1: An example of the distribution of zeros of IZw in the 
complex frequency plane. The vertical axis is In \1ZZ |, so the 
tall spikes represent poles of the scattering matrix. For this 
plot, we have taken (3£ 3/2 = 1.00, ml = 3.34, Hi = 1.37 and 
id = 1.00. 



by cj)*^ and integrate over y. This yields 



dy{WJ 2 + [u 



u 2 H 2 ] 



,| 2 }, (32) 



where we have integrated by parts and made use of the 
fact that (j)^ vanishes strongly as y —> 00 if Imw < 0. 
The lefthand side can be expanded using (|21b[) and the 
boundary condition (|25[) . Then, taking the imaginary 
part of ([32]) . we get 



= Im (w 2 ) 



2(H£f 



dy\cj>u 



(33) 



Since the quantity inside the square brackets is manifestly 
nonzero, we are led to the conclusion Im(o; 2 ) = 0. But, 
we also have that Imw ^ 0, from which it follows that 
Recj = 0. That is, any resonant frequencies u> satisfying 
1Z U = in the lower half of the complex plane must lie 
on the imaginary axis. 



2. Classification of resonances 



where T = Imw. 

We can use this asymptotic behaviour to derive a useful 
result for resonances with Imw < 0. We multiply (|21a|) 



We now turn to the classification of resonant modes 
with complex frequency. One important way of dis- 
criminating between different classes of resonance is 
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the asymptotic behaviour articulated in Eq. pip . An- 
other categorization scheme involves assigning a "Kaluza- 
Klein" mass p to each resonance, which is implicitly de- 
fined by noting that 0^ satisfies the Klein-Gordon equa- 
tion on the brane 



Imw 



-p 2 0« = {d 2 + 3Hd T + k 2 e 
where the field mass is 



- 2HT )C 



p 2 = H 2 (u 



(34) 



(35) 



Note that for Im lu < 0, p will be purely real or imaginary. 
On the other hand, if Imcj > the mass will be complex 
in general. Finally, we can also classify resonances based 
on whether or not they are normalisablc under the Klein- 
Gordon inner product: 



1 



t a {<fd a i> - #>a<n, (36) 



where S T is the spacelike 4-surface of constant r and t a 
is its future directed normal vector. We apply a box- 
normalization of volume V3 for the spatial 3-manifold. 
Consider the norm of <f> R , 



(0S,0S)kg cx / dy 



AH. I 2 



(sinh Hy) 3 



cx / dy\4%\ 2 . (37) 



in, 



Vh 



Hence, </>2 is normalisable if 02 i s square integrable over 
V- 

Armed with these definitions, we can delineate states 
into three categories: 




FIG. 2: Different classes of possible resonances, 
indicate poles of the scattering matrix. 



The circles 



Tachyonic bound states: These have Imw < —3/2 
and Re a; = 0. They are normalisable and diverge 
at late times. In this case, the Kaluza-Klein mass 
is imaginary p 2 < 0, which is why we call these 
modes "tachyonic bound states" . 

We summarize this classification scheme in Fig. [2j We 
have indicated the "mass gap" between p — and p — 
3H/2 familiar from the analysis of uncoupled bulk fields 
about de Sitter branes. 



Scattering states: These are resonances with Imcj > 
0. For these modes, both 02 and 02 vanish at 
late times. Conversely, both states also diverge for 
large y, and the Klein-Gordon norm (02,02)kg di- 
verges. Hence, these are not conventional bound 
states. This spatial divergence does not mean that 
scattering states are unphysical, since they are well 
behaved if the limits r — > 00 and y — > 00 are taken 
simultaneously. Scattering states are akin to quasi- 
normal modes in black hole perturbation theory in 
that they are never observed at spatial infinity, only 
future null infinity. The Kaluza-Klein mass of scat- 
tering states is generally complex. 

Overdamped bound states: These have —3/2 < 
Imw < which implies Rew = 0. At late times, 
the canonical field 02 blows-up while the physi- 
cal field 02 is exponentially damped. These modes 
are normalisable under the Klein-Gordon product 
since (02,0S)kg 1S finite. The Kaluza-Klein mass 
of these modes is real since p 2 > 0. Because 02 
vanishes in the infinite future and is normalisable, 
we call these modes "overdamped bound states" . 



3. Critical coupling and the massless mode 

For any choice of /i, to, Hi > 1, it is always possible 
to choose the coupling parameter (3 = /3 cr it € K such 
that there is one resonance with frequency to = — 3i/2 
and Kaluza-Klein mass p = 0. To see this, we substitute 
lu = —3i/2 in IZ^ = and solve for /3 cr i t : 



/?crit — 



/V2(q-2) 

£1/2 



"(7) 



?f(7) 



1/2 



(38) 



Note that a fundamental property of the associated Leg- 
endre function of the second kind is that [H, 8.783] 

e~" a <X(7) > for all a > -§, a > -|, 7 > 1. (39) 

From this it follows that /3 cr it is real for all possible pa- 
rameter choices; i.e., it is always possible to select the 
coupling such that at least one resonance is massless. 
Note also that (3 CI it = when either to or /i is zero. 

The qualitative behaviour of this massless mode at late 
times is the same a massless scalar propagating in 4- 
dimensional de Sitter: Namely, the physical field ampli- 
tude 02 is frozen (conserved) on superhorizon scales. 
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4- Bound state resonances 

In Fig. [H we track the resonant roots of 1Z W through 
the complex plane as a function of the coupling f3 / /3 cr it 
for particular values of the masses. Here and in Sec. IV CI 
below, the roots are found numerically using Newton's 
method. We see that there are bound states for all cou- 
plings, but there is only a tachyonic mode for super- 
critical (/3 > /3 cr it) coupling. Furthermore, in the small 
coupling regime, we see two overdamped bound states. 
The behaviour of the resonant frequencies for a different 
choice of parameters is shown in Fig. 2] In this scenario, 
no bound state even exists for small coupling. In both 
cases, we see that the fundamental mode goes from over- 
damped to tachyonic as the coupling is increased through 
the critical value. There is no mode mixing in general at 
this coupling, although at other couplings it is possible 
for modes to merge or bifurcate (for example, around 
/?~0.6x /3 cr ;t in Fig. [3}. Our purpose in this subsection 
is to understand some of this bound state behaviour in 
terms of the general properties of 1Z U . 

Massive bulk and brane fields We are only interested 
in the bound states, so we set 



i(A + |), A> 



(40) 



If the Kaluza-Klein mass of a resonance corresponding to 
a given value of A is p\ , then we have 



Pi 



-H 2 X(X + 3) 



{positive, A e [— |, 0), 
negative, A e (0, oo); 



(41) 



i.e., any bound state with A > is a tachyon. 

Assuming that a > 2 and /i > 0, and mindful of (139 
we can re- write the resonance condition as 



(P/Pcrit) 2 = ^(A), 



(42) 



where 



Z(X) = 



1 



(a -2) Q°( 7 ) 



i-4 



(X + a + 1) 



Qxil) 
Qa+i(7) 



7(A + 3) 



(43) 



Empirically, we find that for all a > 2, fi > and 7 > 1, 
the continuously differentiable function Z(X) € K has the 
following properties: 

• Over the interval A E [0,oo), Z(X) monotonically 
increases from 1 to 00. Therefore, there always 
exists one and only one tachyonic bound state for 
supercritical coupling (3 > /3 C rit- 

• In the interval A 6 [—3/2, 0], Z{X) has at most one 
local minimum and no local maxima. Therefore, 
for any choice of parameters we can have at most 
two bound state resonances. 



• The absolute minimum of Z(X) for A € [—3/2, 00) is 
always less than unity, but not necessarily less than 
zero. Therefore, it is possible to select parameters 
such that there are no bound states at all. Such 
situations will always have subcritical coupling (3 < 

/3crit • 

We reiterate that these conclusions are only valid with 
both fields have finite mass. 

From these comments, we can deduce that a neces- 
sary (but not sufficient) condition for the existence of 
two bound states is 



Z(-3/2) > and Z'(-3/2) < 0. 



(44) 



By direct sampling of the (fj,/H,m£, Hi) parameter 
space, we have numerically determined that these con- 
ditions hold if and only if 

fj,/H < 3/2 and m 2 /H 2 < m 2 hl . cs (H i) . (45) 



Here, the dimensionless threshold mass rh 2 hTes (H£) is ob- 
tained numerically and shown in Fig. [5] We will see in 
Sec. IIII C II below, that the second inequality in (|45|) is 
a necessary and sufficient condition for there to exist a 
bound state of cj) at zero coupling. 

Massless bulk field The situation is somewhat differ- 
ent when either m or /i is set to zero. For the sake of 
brevity, we only mention the case of the massless bulk 
field here. When a = 2, the condition for bound states 



f3 2 /H 3 =y(X), 



where, 



y(x) = 



2p 2 > 2 A -/i 2 ) 

I/ 4 



(46) 



(47) 



The function y(X) € R is continuously differentiable for 
A > —3/2, and has the following properties: 

• y(X) monotonically increases from to 00 for A 6 
[0, 00). Therefore, there always exists one and only 
one tachyonic bound state for any nonzero coupling 
(3>0. 



• The sign of y at A = -3/2 satisfies sgn [3>(-3/2)] = 
sgn [%H — fi]. Therefore, in addition to the ubiqui- 
tous tachyonic mode, one can have an overdamped 
bound state if and only if the mass of the brane 
field lies in the mass gap /i < §7?. 

• y(X) has at most one local extremum for A > —3/2. 
If this extremum exists, it must be a minimum in- 
side [—3 /2, 0] . Hence, we can either have one tachy- 
onic and one overdamped resonance, or one tachy- 
onic resonance only. 
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(a) A few of the complex roots of Ti^j as (b) Variation with /3//3 cr it in the imaginary parts of the 

functions of the coupling parameter f3 resonant frequencies shown in (a) 



FIG. 3: The least-damped resonant frequencies of the system when fi£ = 1.40, ml = 2.37 and HI = 1.25. This particular 
example illustrates the generic feature of our model: a tachyonic resonance only exists for supercritical (f3 > /3 cr it) coupling. 
Note that all resonances with Imu < have Kecu = 0, as expected. Also note that the two /? = resonances uo — ±0.998i 
are modes for which the bulk field vanishes cj> u = (i.e. these are pure q modes), while the other modes have a vanishing 
brane mode q^ = 0. The "fundamental mode" has been identified as the resonance with the smallest imaginary part, which is 
expected to dominate the late time behaviour of numeric simulations (cf. Sec. IV C|) . 




PI flciit FIG. 5: The threshold mass as obtained from numeric cal- 

culations. The condition m 2 /H 2 < rh 2 hTCE (H£) is necessary 
FIG. 4: Variation in the imaginary parts of the least-damped for the existenc e of two bounds states when the coupling is 
resonant frequencies with coupling when fd = 0.40, mi = 2.1 nonzero. In Sec. MIC II we also see that this condition is nec- 
and Hi = 0.20. In contrast to Fig. O there is no bound essary and sufficient for the existence of a bound state of the 
state at small coupling. For /3//3 crit > 0.65, there is a single bulk field at zero coupling, 
bound state that is overdamped or tachyonic for subcritical 
or supercritical coupling, respectively. 

1. Zero coupling 

We first consider the zero coupling limit of the model 
C. Limiting cases /3 = 0. It is useful to define 

ftW=(a-iu,-l/2)Q? +1/2 ( 7 )+ 

In this subsection, we consider the resonance condition ^ (f w — 3 /2) Qf u _ t / 2 (7), (48a) 

TZaj = in various limits and hence recover several results , , „ „ „ 

from the literature. ^ = t* l H ~ " ~ 9 / 4 - ( 48b ) 
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In terms of these functions, we find that 1Z U = implies 

7#> = or 7#> = 0, [/3 = 0]. (49) 

Going back to the analysis of Sec. lIII Al we can show that 
if either the first or second condition holds, then <7^(t) = 
or 0S( T i2/) = 0, respectively. In other words, each 
condition represents the resonances of one field decoupled 
from the other. The brane field resonances are trivial 

,2. 



= ±vV/ff 2 -9/4 



(50) 



p =n 

that is, these modes oscillate with the natural frequency 
of the brane oscillator in isolation and have a Kaluza- 
Klcin mass equal to p. 

Using an identity [23], 8.739] and recurrence relations 

(23l 8.73], the bulk resonance condition izffl — can be 
re-written as 



= C(« 
= (iu - 



a ■ 



1\ p—iu) 
2 > a+1/2 



(C) + (i«-a + |)P-^ /2 (C) 

(C)-C(« + 2)P-^ /2 (C), (51) 

where £ = j/(H£). These match the various bound state 
criteria obtained for a massive scalar field about a de 
Sitter brane in Refs. [13, [H, [29| (up to trivial differences 
in notation). 

As in Sec. MI B 41 above, we can substitute in ui = 
+ |) into the bulk resonance condition 7?if^ 
to derive a criterion for bound state resonances of 
zero coupling. This is 

0a +2 (7) 7 A 



= X{\) 



QS +1 ( 7 ) \ + 2-a 



A > 



= 
d> at 



(52) 



One can confirm that X(X) £ 1 is monotonically in- 
creasing for A > —3/2 and positive definite for A > 0. 
Therefore, a de-coupled massive bulk field can exhibit at 
most one bound state, and this must be an overdamped 
resonance. Furthermore, a necessary and sufficient con- 
dition for the existence of such a bound state is that 
X{— 3/2) < 0, which we have numerically determined is 
equivalent to 

m 2 /H 2 < m 2 hlcs (H£), (53) 

where the function m 2 hres (iJ£) is shown in Fig. [5] In 
other words, if the bulk field is too heavy compared with 
the Hubble scale, it cannot be effectively localised near 
the brane. 

Finally, setting the bulk field mass to zero (i.e. a = 2) 
in (|48ap yields that 

-10 31-1/2(7)- (54) 

In this case there is always a critical bound state with 
lu = — 3z/2, which is usually called the zero mode; i.e., 

zero mode to = —3i/2, m — 0, (3 = 0. (55) 

One can confirm that for the zero mode, the resonant 
solution is y-independcnt 



cV> R (r,y) =0. 



(56) 



This fact will be useful in the interpretation of the effec- 
tive 4-dimensional theory in Sec. HVl 



2. Perturbed zero mode 

Since the coefficients in front the Legendre functions 
in (|29|) vanish when m = and w = — 3i/2, it is possible 
to linearise the resonance condition about the zero mode 
solution (|55p . In particular, we can take 

(57) 



|w + |i|<l, m£<l, 10 2 < min(^ 2 ,i? 2 



expand 1Z U to linear order in small quantities, and solve 
for u. Substituting this solution into (|35[) . we obtain 

SpLo - m 2 g(H£) - £~ x ^ 2 (i 2 f{H£) 

~m 2 g(H£)(l-p 2 /p 2 rit ), (58) 

where Sp zcm is the small Kaluza-Klein mass that the zero 
mode acquires when either m or (3 is nonzero. Here, we 
have defined 



1 - x 2 In 



1 + Vx 2 + 1 



^ 3 



(59a) 



(59b) 



and made use of the limiting values of the critical cou- 
pling ((38]) for m£ < I : 



Petit 



£m pt 



2..2 9(H£) 



f(H£) 




\m 2 p 2 £, 
2 m 2 p 2 H-\ 



H£<€1, 



(60) 



as well as explicit representations of the Legendre func- 
tions with integral indices [23], 8.827]. We see that the 
mass of the bulk field tends to increase the squared mass 
of the zero mode into the overdamped region, while the 
coupling tends to decrease it into the tachyonic regime. 
We also see that supercritical coupling (3 > /3 cr j t explic- 
itly leads to a tachyonic mode 5p 2 cro < 0. For reference, 
we note the limiting behaviour of Sp z 
on extreme values: 



when H£ takes 



dp: 



zero 



\m 2 
|m 2 



£r L p- 2 f3 \ H£<1, 
lp- 2 Hf3 2 , H£^>1. 



(61) 



Note that for the perturbed zero mode to be a bound 
state, we need 

Sp ZCIO < 3H/2 => m/H < V3/2, [HI < 1]. (62) 

This matches the small H£ behaviour of the bound state 
condition m 2 /H 2 < w 2 hrcs , where m 2 hrcs is shown in 
Fig. There is no conclusion to be drawn in the other ex- 
treme, because the conditions (|57[) imply 5p zeTO /H <C 1 
for H£ 3> 1; i.e., the perturbed zero mode is always a 
bound state. 

We note that the expression (|58j) for Sp ZCIO closely re- 
sembles that obtained by Langlois and Sasaki [2^], who 
studied a massive bulk scalar with a brane-localized self- 
coupling. Their model was governed by the action 



S hS = - (g ab d a ct>d b( t> + m 2 <t> 2 )- 



(63) 



M 



8Mb 
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where x is a dimensionless constant parameter. Their 
results for the perturbative zero mode mass match ours 
if we make the identification 



(64) 



Hence, under the conditions (|57|) . our system behaves 
like a massive bulk scalar field with a constant brane 
self-coupling. 



3. Low energy 

The low energy limit of the model is defined by 

#£<1, 7-1, [£ finite]. (65) 

Consider the temporal eigenvalue equation (fT?)) in this 
regime: 

(<9 2 + fc 2 )TL,(r) = -Wp hys T w (r), w phys = u>H, (66) 

where we have identified the dimensionful physical fre- 
quency Wphys- When we take the limit Hi — > 0, we want 
to ensure that the physical frequency is finite. Hence, we 
define p 2 /i 2 to be the finite eigenvalue of (d 2 + k 2 ): 

p = luH£ = £iv D h vs , lim \p\ < oo, lim \u>\ — oo. (67) 

1 y H£->0 Hi-iQ 



In terms of p, the solution of (|66|) is 

T w (t) = e mT '\ w = ^/p 2 + k 2 i 2 . (68) 

To find the resonant modes in this case, our strategy will 
to be to take the limit of IZ^ — when 7 — > 1, Hi <C 1, 
I a; I ^> 1, and w real. To avoid branch cut ambiguities, we 
put p — > p — i0 + with p € M. The resulting expression 
can then be continued into the p complex plane. 

Consider the following integral representation of the 
associated Legendre function [23|, 8.711(4)]: 



E(7) 



e i7rQ T(er + 1) 



dt cosh at 



T(a-a + l)J (7 + ^ 7 2 -lcosht) ff + 1 
[7 > 1, Rc(cr-a) > -1, and -a ^ 1,2,...]. (69) 

Now, if we assume 7 — > 1 and a = iui ± 1/2 — i(/3 — 
zO+Xf/TT 1 in the integrand of {69|), we obtain 



e i7r T(w + 3/2) 



r(iw- ^ + 3/2) 

/•OO 

/ dt exp [- i (p - i0 + ) cosh t] cosh i4 +£>(#£) , 
Jo 



= -ii7re--"/ 2 H, ( , 2) (p-iO+) 



(70) 



where the integral [23|, 3.457(4)] has been evaluated as- 
suming the principal branch of the Hankel functions 
H^' 2) (z) [H 8.476(7)] and fundamental definitions [H, 



8.407]. Using ([6"7 ) -([?0 |) . a recurrence relation [23j, 
8.734(4)], and expanding to leading order in Hi, we ob- 
tain the resonant condition for low energies: 







■l^ + (a-2)(p 2 -p 2 i 2 )] H( 2 >(p-*0+) 



p{p Z ~ ^^Up M + ) + O(Hi). (71) 



This matches the result obtained in Koyama et al. [l6| if 
one takes into account different conventions for the field 
time dependence and an implicit choice of integration 
contour for the inverse Fourier transform. Note that in 
this limit, the Kaluza-Klein mass (|35[) of a given reso- 
nance is 



9H 2 m ; j? 
4 t 2 ■ 



(72) 



Hence, the critical massless mode corresponds to p — 
0. As before, this mode is a resonance if (3 equals some 
critical value 



Aait =^" V V2(a-2). 



(73) 



As discussed in detail in Ref. [16| . bound states in this 
limit correspond to resonances with p on the negative 
imaginary axis. One can only have a bound state for 
supercritical coupling, and that unique bound state is 
necessarily tachyonic. 

Eq. ([71]) can be linearised about the resonance {55]), 
resulting in the mass perturbation of the zero mode: 



Sp 



2 

zero 



1 2 

— m 

2 



p z i 



1 2 
— m 

2 



1 



8 2 



(74) 



where the last approximation comes from expanding the 
critical coupling about m = 0. One can confirm that this 
reproduces the Hi <C 1 limit of Eq. (f58f . 

We can also take the limit of ([71]) in the case of a flat 
bulk 



I — » 00 , p = pi ~ > 00 , mi ~ a 
by making use of the following fact [Til ]: 



(p, Q ) 



\J m 2 i 2 — p 2 , 



(75) 



(76) 



HL 2) (P) ' 00 

We obtain 

— P 2 + 2\/m 2 ~~p 2 (p 2 — p 2 ), p = p/i, (77) 

which matches the resonance condition found by George 
[I5I ] for the case of a Minkowski brane in a Minkowski 
bulk (with an appropriate change of notation p 2 — > p 2 + 
m 2 and choice of branch y— p 2 = — ip). 

^. ffig/i energy 



The high-energy limit is 

7~JI^>1. 



(78) 
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Unlike the low energy limit above, we make no assump- 
tions about the magnitude of ui. Using the asymptotic 
expansion for Q" LU ± 1 / 2 (l) f° r large 7 [23f, 8.776(2)], and 
working to lowest order in (Hi)^ 1 , we find that the res- 
onant condition 1Z W = reduces to: 



fP_ 1 
H 3 ~ iuj + 1 



2u 2 



9_]J?_ 
4 H' 2 

\p/H and mi non-infinite]. (79) 



m 2 \ 



Note that the resonance condition is independent of £, 
which is expected because the characteristic curvature 
scale on the brane is much smaller than the char- 
acteristic curvature scale of the bulk £; i.e., the brane 
"sees" the bulk as flat. 

Eq. (fT9")l implies that there are four resonant frequen- 
cies that are the solutions of a quartic polynomial. The 
critical coupling is again obtained by setting uj = —3i/2: 



(80) 



which matches the Hi — > 00 limit of Eq. (f38|) . As 
in Sec. Iff! C 21 we can linearise the resonance condition 
about the zero mode (|55|) , resulting in the perturbation 
of its Kaluza-Klein mass: 



5p; 



3 o 3// n 3 y ( (3 2 
— ~ 5 m " = 5 m I 1 " ^ 



(81) 



This matches the large iff limit of ([55)1. as expected. 

Note that the high energy limit (|78|) does not distin- 
guish between H approaching infinity or i approaching 
infinity. Hence, it is possible to expand (|79| in the case 
of a flat stationary brane; i.e., H — > and i 00. This 
yields 



(3 2 = ip~ 1 (p 2 - fj?)(m 2 -2p 2 ), P = ujH. 



(82) 



Note that this result implicitly assumes the m — ► limit, 
since (|79[) is only valid when mi is non-infinite. It is 
easy to confirm that this matches the flat space resonance 
condition obtained in the low energy regime (|77[) when 
m — > 0. 



IV. 4-DIMENSIONAL EFFECTIVE THEORY 
AND THE BOUND STATE MASS MATRIX 

In this section, we show how many of the qualitative 
and quantitative properties of the model derived from 
the spectral analysis in Sec. IIIII are encapsulated by an 
effective 4-dimensional action. Here, <j> and q refer to 
the fields in real space, rather than the spatial Fourier 
transforms. 

To derive this action, let us write 

0(r,x,y) = pi(<r,x)f(y), g(r,x) = ^ 2 (r,x). (83) 



We substitute these definitions into the full 5-dimensional 
action @ and explicitly perform the integrations over y. 
We obtain the following 4-dimensional effective action: 



Scff 



2 

E 

i=l 



1 2\ 



2 (4) /3^i^2 



(84) 

Here, Greek indices are raised and lowered with the brane 
metric h a p. We have defined the effective masses and 
coupling as 



2m' 



dyN 5 (y)e(y) 



2 / dyN 3 (y)[Z'(y)} 2 , 



ml 



where N(y) — Hi{smhHy) 1 
assumed 



<4) /3 = fc{y h ), 
In deriving 



(85) 
we have 



1 

dyN 3 (y)e(y) = 



This assumption is equivalent to demanding that the 
original bulk field is normalisable under the Klein- 
Gordon inner product (f36|) . Hence, our effective action 
only gives a valid description of the bound state degrees 
of freedom in the model. The full action ([9]) is required 
if we want to include continuum modes (i.e., scattering 
states) as well. 

The 4-dimensional action (|5^)) represents a pair of cou- 
pled massive fields confined to the brane. We can seek 
mass eigenstate solutions (i.e., eigenfunctions of V a V„) 
as follows: 



n\ < 4 >/3 
»/3 ml 



M 



Mxf> = p 2 Tp 

\ H>2 



(87a) 
(87b) 



This matrix equation has a non-trivial solution if p 2 is 
one of the eigenvalues of the mass matrix M, which are 



P± = 



1± 



4(< 4 >/3 2 - m\m 2 ) 



(mj 



(88) 



These are the masses of the two normal modes of the 
system, which should be identified with the Kaluza- 
Klein masses of bound state resonances in the full 5- 
dimensional theory. 

The fact that we can have at most two mass eigen- 
states in the 4-dimensional effective theory is consistent 
with the results of Sec. MI B 4( where we saw that there 
can be at most two bound state resonances in the full 
5-dimensional treatment. Furthermore, we note that one 
of these modes becomes tachyonic (i.e., pi < 0) if 



(4)o2 



■f3 z > w f3 2 



(89) 



This again matches our expectations: We have one and 
only one tachyonic mode if the coupling is sufficiently 
large. 



11 



It is important to note that not all mass eigenstates 
in the effective theory spectrum correspond to bound 
states of the full 5-dimensional model. Recall that bound 
states must have Kaluza-Klein mass within the mass gap, 
so any mass eigenstates with p\ > 9H/4 are not true 
resonances. We have that p\ is bounded from below 
by max(mj,m|), while p_ is bounded from above by 
min(r7i 2 , m 2 .). Thus, a necessary condition for there to 
be two bound state solutions in the effective theory is 

mi < 3H/2, m 2 = p< 3H/2. (90) 

Qualitatively, this matches the results of Sec. IIII B 41 
where we saw that a necessary condition for the exis- 
tence of two bound states was that the field masses were 
not too heavy. Unfortunately, one cannot translate these 
inequalities into more stringent constraints on m without 
knowing £(y). We should briefly mention the Minkowski 
brane case: When H = there is no mass gap and only 
tachyonic resonances with p 2 < correspond to bound 
states. Hence, in the effective theory there is only a sin- 
gle tachyonic bound state for supercritical coupling, in 
agreement with the conclusions of Ref . [l(| • 

We can quantitatively compare our effective theory 
to the results previously obtained for the perturbed 
zero mode in Sec. IIII C 21 We assume that the extra- 
dimensional profile is the purely "zero mode" solution, 
as given by Eq. (|55|) : 

£(y) = constant. (91) 

Under this assumption, the various integrals in (|85|) and 
(1551) can be reduced to elementary functions. We find 

e = r\f(H£), ml = m 2 g(H£), 

=r 1 p 2 f(H£), (92) 

where / and g are defined in Eqs. (|59|) . To relate these 
back to the perturbed zero mode, recall that our re- 
sults for 8p 2 cvo were derived for the case when m and 
(3 were small. This suggests that we expand the effec- 
tive theory resonant masses (|88[) under the conditions 
m| 3> max(m 2 , <4) /3), which yields 

p\ ~ m\ + w f3 2 /ml p 2 _ ~ ml - {i) p 2 /m 2 2 . (93) 

Substituting (|92p and m 2 — p into the above expression 
for p 2 _, we recover the the 5-dimensional prediction for 
<!>Pzero Finally, we note that the effective theory 

prediction of the critical coupling (|89|) in this regime is 

< 4) /3 c 2 rit = p 2 m 2 g(H£) ~ r 1 f(H£)f3 2 clit , (94) 

where we have made use of (|60| . The relationship be- 
tween <4) /3 2 rit and /3 2 rit is entirely consistent with what 
we expect from (|9"2"1) . 

To summarize this section, we have shown that the 
bound state sector of our model can be described by the 
effective action of a pair of coupled scalars propagating 
on a de Sitter background. The main complication is 



that the effective 4-dimensional masses and coupling de- 
pend on higher-dimensional details, so moving from the 
5-dimensional to 4-dimensional pictures requires the non- 
trivial calculations performed in the rest of the paper. 

V. NUMERIC TIME DOMAIN ANALYSIS 

We now turn to the solution of the coupled equations 
([TO)) subject to the boundary condition (jll| in the time 
domain. As in Sec. IIII! we will work with the spatial 
Fourier transforms (|12p of </> and q and omit the k sub- 
script. 

A. Dimensionless wave equations in Poincare 
coordinates 

In our numeric integrations, we elect to use dimension- 
less versions of the Poincare coordinates t and z (J2J), as 
well as the conformal time r\ ([3]). In particular, we define 

t = t/z*, z = z/z*, f) = r)/z*, 

B t = d/di, d z = d/dz, d n = d/dfj, (95) 

where is a length scale we specify later. We rescale 
the bulk field as follows 

4>(t,z) =£ 1/2 z- 3 / 2 (t)(t,z), (96) 

and write the brane trajectory as 

tb = -ffj, % = —Hlfj, dfj 2 = dt^ — dz 2 . (97) 

Then, the wave equations (|10p and boundary condition 
(|TT|) reduce to 

= [d 2 - d 2 + V(z)]ip, (98a) 
- [d 2 + A(z h )d v + B(z h )]q + C(z h )^ b , (98b) 
= (cW)b - D(z h )i/j h - E{z h )q - F(z h )d v q, (98c) 

where the potential is 

V(z) = k 2 z 2 + ±(15 + 4m¥ 2 )z~ 2 , (98d) 

the coefficients are 

A(z h ) = 2H£z- 1 , D(z h ) = -Ijz-\ 

B(z h ) = k 2 ? 2 + pH 2 ?^ 2 , E(z h ) = \(3P/ 2 z- 5/2 , 
C(z h ) = 0£ 3 /% 1/2 , F{z h ) = 0, (98e) 

and the (flat-space) normal derivative operator is 

d n = -Hld t + 1 d z . (98f) 

Note that even though the F coefficient is zero for the 
problem at hand, we retain it to keep the algorithm de- 
veloped in the next section reasonably general. 
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141 



bulk cell 
(diamond) 





boundary cell 
(triangle) 



FIG. 6: The spacetime domain Q over which we seek a nu- 
meric solution of (I98|l . Superimposed on Q is a (particularly 
coarse) example of the computational grid we use to discretize 
the problem. Details of the two types of cellular geometries 
are shown in the upper right corner. 



B. The algorithm 

Previously, Seahra [IH has developed a numeric algo- 
rithm to solve for ip in equations similar to (|98|) when 
there is no brane-bulk coupling (/? = 0) and the brane 
trajectory is arbitrary. Here, we generalize that proce- 
dure to the case when there is a dynamical degree of 
freedom on the brane. 

We seek the numeric solution for ip throughout a finite 
region of the (t, z) plane O depicted in Fig. [6j as well as 
the solution for q on a finite segment 9f2b of the brane 
trajectory. To have a well-posed Cauchy evolution for 
the fields, we must specify initial data for ip on the past 
null boundary of the computational domain dfl~, as well 
as the value of q and d v q at V = <90b n dtt~ . 

As in Ref. [18(, our strategy will be to partition ft 
into a large number of small cells to construct a discrete 
approximation to the fields. In Fig. [51 we see that the 
cells bounding the brane take the shape of triangles with 
one timelike and two null sides, while the cells interior to 
Q are diamonds with equal null sides. By integrating the 
bulk wave equation ()98a[) across a diamond cell, one can 
find out the value of ip at the northern node of the cell 
(i.e., the uppermost corner) in terms of its value at the 



other three nodes (cf. 18]): 

1>n = -i's + + i'e){l - \5 2 t V s ) + O(Sf), (99) 

where St = t n — t s is the difference between t evaluated 
at the top and bottom of the cell. The subscripts n, s, e 
and w indicate the value of a quantity at the northern, 
southern, eastern and western nodes, respectively. 

The evolution across a triangular cell is more compli- 
cated due to the presence of the boundary degree of free- 
dom q. It is useful to re-write the second-order ODE 
(|98b[) as a system of two first-order equations 

dr,q = p, = d, lP + A(z h )p + B(z h )q + C(z h )ip h . (100) 

These can be discrctized across a given triangular cell 
using the modified Euler method: 

q n = q s + 5* n (p» + P.) + G(S*), (101a) 
Pn = Ps~ \S 71 (A n p n + A s p s + B n q n + B s q s + 

+C n lP n + C s lP s ) + 0(S*). (101b) 



Also, by integrating (|98ap over the cell and making use 
of (|98cP we obtain 

6S v (E n q n + E s q s + F n p n + F B p s ) 



1pn = - 



12 + 6S v D n + S v„ 
12 + 6S V D S +S U 5 V V S 



12 + 65 n D n + S u S v V n 
2A-S u S v V e 
12 + Q5 n D n + S u S v V n 



4> s 



(102) 



Here. 



St - 8 Z , 
St 



S v = (Sl-S 2 z )V\ 



(103) 



By solving the linear system p01[) and (|102[) . we can 
obtain the future field values (q n ,p n , ip n ) in terms of the 
past field values (q s ,p s ,ip s ,ip e ). 

Having written down these evolution formulae, the ac- 
tual algorithm for obtaining tp and q is very similar to the 
one used in Ref. [l8j. Consider the first row of cells in 
Fig. [51 If we know the value of -0 along dil~ as well as the 
value of (q,p) at V, then application of the triangle and 
diamond evolution formulae in sequence (i.e. southwest 
to northeast) allows us to obtain (q,p^ijj) on the future 
boundary of the first row. One can then find the values 
of ip on the nodes of the past boundary of the second 
row via interpolation. By repeating this procedure for 
rows 2 through AT, one obtains the solution for the fields 
throughout ft. As in Ref. (l8j . we expect our numeric 
solutions for ip and q to converge like Sf to the actual 
field values. 
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C. Results 



1. Minkowski brane 



We now apply our algorithm to the problem introduced 
in Sec.|TTJ First, let us consider the simpler case of a brane 
stationary with respect to the Poincare frame: 



HI = 0, z h = I, 



dsl 



-dT l + dx 2 



(104) 



i.e., a brane with Minkowski geometry. In this case, we 
select 



z* = t 



1. 



(105) 



As discussed in Sec. IIII C 3[ the time dependence of any 
given resonant mode in this case is given by 



(106) 



where p is a solution of the low energy resonance condi- 
tion Ijnj with H£ = 0. 

As discussed in Sec. IIII B 11 we expect the late time 
behaviour of <f> and q to be dominated by resonant mode 
solutions at late times, no matter what initial data is 
selected on dQ~ . Furthermore, if we wait long enough 
the resonant mode with the smallest imaginary part — 
which we call the fundamental mode — will dominate 
the other modes. Hence for any choice of initial data, we 
expect that 



q ► Re (C„e v 

OO 



Re(C 



), (107) 



where C q and are complex numbers, and vjq is the 
frequency of the fundamental mode. 

In Fig. we show the results of our numeric simula- 
tions for the following choice of initial data: 



q(V) = l, p{V)=Q = (j>{dQ.-). 



(108) 



We find excellent qualitative and quantitative agree- 
ment between the simulation results and the expecta- 
tions for the fundamental frequency from the spectral 
analysis. In particular, we get exponentially diverging 
fields (i.e. tachyonic resonances) for large coupling. We 
further test the agreement between simulations and the 
resonance condition (|71j) in Fig. [H 



2. De Sitter 



For a de Sitter brane, we select to be the position of 
the brane when the mode under consideration exits the 
Hubble horizon, which implies 



kz„ = He. 



(109) 



Along the same lines as our discussion of the Minkowski 
brane above, we expect that the late time behaviour of 



w pied = 1.13 + 0.00363* 
w sim = 1.13 + 0.00362? 
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FIG. 7: Simulation results for a Minkowski brane with \xl = 
1.23, ml = 2.37 and ki = 0.340 and the initial data ([1081) . 
The quoted values of ro plC d are predictions for the fundamen- 
tal frequency from the low energy resonance condition (I71|l . 
while tj7 aim is the late time frequency obtained by fitting the 
template (|107p to the numeric waveforms. This plot exhibits 
the three possible classes of asymptotic field behavior when 
H£ = 0. 



the fields in a numeric simulation to be 



Re 



6b > Re 

OO 



^ e (iu -3/2)Hr 
q (io) D -3/2)Hr 



(110a) 
(110b) 



regardless of the particular choice of initial data. Here, 
the fundamental frequency uj$ is the solution of (|29[) with 
the smallest imaginary part. 

In Fig. we show results for some of our simulations. 
The particular case considered is one in which the funda- 
mental mode is a scattering state for /3//3 C rit ^ 0.65, and 
an overdamped bound state for 0.65 < (3/f3 cr it < 1- In 
Fig. 1101 we quantitatively compare the value of the fun- 
damental frequency obtained from the spectral analysis 
and the value obtained from simulations in a case where 
there is an overdamped bound state for all subcritical 
couplings. We find acceptable agreement between the 
values of luq obtained in both approaches. 
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FIG. 8: Comparison between spectral analysis predictions for 
the fundamental frequency and simulation results when fil = 
1.23, mi = 2.37 and hi = 0.340 (Minkowski brane). 
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FIG. 9: Simulation results for a de Sitter brane with /j£ = 
0.40, mi = 2.1 and Hi = 0.20 and the initial data (fT08l> . Here, 
Ti, is the proper time coordinate on the brane when the mode 
exits the Hubble horizon. Individual curves are labeled by the 
coupling parameter and the classification of the fundamental 
mode solution of (|29[1 for these parameters — the four possible 
varieties of late time behaviour when HI > are all seen in 
this plot. 
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FIG. 10: Comparison between spectral analysis predictions 
for the fundamental frequency from the full resonance con- 
dition TZui ~ and simulation results when fj£ = 1.00, 
mi = 2.24 and Hi = 1.12 (de Sitter brane). For this choice of 
parameters, the fundamental resonance is never a scattering 
state. 



VI. CONCLUSIONS 

In this paper we have studied resonant modes of a mas- 
sive bulk field in five-dimensional anti-de Sitter spacetime 
linearly coupled to a massive field on a four-dimensional 
boundary de Sitter brane. Using spectral scattering the- 
ory in 5 dimensions, we investigated the resonant exci- 
tations of the system, identifying scattering states which 
propagate to future null infinity and normalisable bound 
states. By dimensionally reducing the full action, we have 
also derived an effective 4-dimensional theory that gives 
a valid description of the bound state resonances. The 
effective theory is essentially that of a pair of coupled 
massive brane-confined scalars. 

From the 5-dimensional treatment, we find that the 
bound state sector of the spectrum is characterised by 
zero, one, or two resonances. These can be divided into 
two categories: overdamped and tachyonic modes which 
decay or diverge exponentially at late times, respectively. 
We find that there is a critical coupling (3 clit , above which 
there always exists a unique tachyonic mode; conversely, 
no tachyonic bound state exists for subcritical coupling. 
The critical coupling is given analytically by Eq. {38}, 
and we find that f3 cr n = if either field is massless. 

All of these bound state properties are mirrored in the 
effective 4-dimensional theory: namely, the effective ac- 
tion predicts the existence of at most two mass eigen- 
states, or normal modes, which correspond to bound 
states of the higher-dimensional theory. Furthermore, 
one of the normal modes is tachyonic when the coupling 
exceeds a critical value, and that critical value is zero if 
either of the effective field masses vanish. 

One of the more interesting features of our model is 
the existence of overdamped bound states, which do not 
appear in the case of a Minkowski brane in an AdS bulk 
[lq . It is well known that for a Minkowski brane, all bulk 
modes with real Kaluza-Klein mass reduce to travelling 
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waves far from the brane and are hence non-normalisable. 
However for a de Sitter brane, it is possible to have nor- 
malisable bulk modes with Kaluza-Klein mass less than 
2>H/2 [l7| . The overdamped bound states found in this 
paper live within this mass gap, which explains why they 
are not present when the brane is "at rest" (H = 0). 
For these stationary branes, one can only have a single 
tachyonic bound state, and that state only exists for su- 
percritical coupling. 

The actual number of bound states for a given choice 
of parameters is strongly influenced by the masses of the 
bulk and brane fields. In particular, a necessary condition 
for the existence of two bound state resonances is that 
both fields are "light" compared to the Hubble scale. In 
this context, the notions of light and heavy refer to the 
behaviour of the individual fields at zero coupling. The 
brane scalar is considered to be light if it does not oscil- 
late on super-horizon scales when f3 = 0. Quantitatively, 
this means that the boundary mass is less than 3H/2. 
On the other hand, a light bulk scalar is one for which 
there exists a normalisable bound state in the absence 
of any brane-localised couplings. We have numerically 
determined that this is equivalent to demanding that the 
bulk mass squared m 2 is less than ~ 4.50 H at low en- 
ergy {HI <C 1) and less than ~ 4.42 H at high energy 
{HI 3> 1). We can qualitatively recover these necessary 
conditions for multiple bound states from the effective 
theory, but we can only obtain the precise bound on m 
quoted here from the full spectral formalism. 

We studied the behaviour of the model's resonances 
in a number of different limits in Sec. IIII C[ and have 
recovered several results previously available in the liter- 
ature. One case worth emphasizing involves the pertur- 
bation of the massless de Sitter zero mode when (3 and 
m take on small values. Under these conditions, we have 
demonstrated that there is complete agreement between 
the predictions of the 5-dimensional spectral formalism 
and the 4-dimcnsional effective theory for the pertur- 
bative zero mode mass, critical coupling, etc. In par- 
ticular, we find that the perturbed zero mode mass is 
given by Sp 2 zcro ~ \m 2 {\ - /3 2 /#? rit ) for HI < 1 and 
8pl eto ~ fm 2 (l - /?7/? c 2 rit ) for HI > 1. From the con- 
dition that the the Kaluza-Klein masses of normalisable 
states must lie within the mass gap, we have confirmed 
that a necessary condition for a bulk bound state at zero 
coupling is m 2 / H 2 < 9/2 for Hi <C 1, in agreement with 



the numerical results of Sec. IIII C 11 

We have tested the predictions of the spectral theory 
against a new numerical algorithm that is able to evolve 
the coupled brane and bulk fields. We are able to con- 
firm our spectral analysis predictions for the fundamen- 
tal frequencies by studying the late-time behaviour for 
both a static (Minkowski) brane and a moving de Sit- 
ter brane in the Poincare coordinates. In particular we 
find good quantitative agreement between the rate of de- 
cay or growth of the boundary field and the imaginary 
part of the frequency of the fundamental mode for both 
subcritical and supercritical coupling. 

We have chosen to work with a moving brane Poincare 
coordinates so that the code can also be used to study 
cosmological branes with arbitrary trajectories in AdS. 
Our numerical code can be used to study the coupled 
evolution any system where the bulk and boundary wave 
equations and boundary condition can be written in form 
given in Eqs. ([98]) . We will thus be able to study cou- 
pled bulk and brane fields in situations where we have 
no analytic results to compare against. For instance 
such a numerical code is required to study the evolution 
of cosmological density perturbations consistently cou- 
pled to bulk metric perturbations during the matter- or 
radiation-dominated eras. 

The simple model of linearly coupled bulk and brane 
fields on a de Sitter brane which we have considered in 
this paper may be a useful toy model for the study of 
inflaton perturbations linearly coupled to scalar metric 
perturbations during inflation in a higher-dimensional 
bulk. However the inflaton-metric interaction involves 
derivative couplings which makes the problem consider- 
ably more challenging. In particular we find that the 
bulk field is coupled to an infinite tower of bulk modes 
[l2l Il3l |30| . We will return to this problem in future 
work. 
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